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The quantization of the electromagnetic field in the presence of material bodies, at zero tempera- 
ture is considered. It is shown that a dielectric does not act as thermal bath for the field and yields 
a non-trivial non-thermal mixed state of the field. The properties of this state and its entropy are 
studied. The dependence of the second Renyi entropy of the field on the distance between dispersive 
objects is shown to decay as R~ 4 for generic bodies. 



INTRODUCTION 

The study of quantum correlations and entanglement 
in many body quantum systems has become a major 
theme in recent years. In particular, many-body sys- 
tems are studied in terms of the reduced density ma- 
trices describing the state of constituents of the system. 
This point of view gives rise to various questions such 
as the question of thermalization: When can a part of a 
quantum system be considered as effectively in a thermal 
state? 

Naturally, thermal density matrices appear in systems 
weakly coupled to a thermal heat bath. Moreover, it has 
recently been argued that a canonical thermal state may 
arise for typical pure states of the system+bath Hilbcrt 
space after tracing out the bath [TJ [5]. However, many 
physical systems of interest are not typical and violate 
some of the assumptions leading to thermalization. In 
particular, the system considered here, that of a boson 
radiation field interacting with a dielectric medium, char- 
acterized by a dielectric function s(uj), has been shown 
[5] to be described by a density matrix which is not ther- 
mal for the typical frequency dependence of the dielectric 
function e(w). 

The dielectric function serves as a convenient way to 
encode the particular long wavelength properties of the 
material [3] . The dielectric susceptibility of the material, 
and it's relations with various response functions serves 
to define an effective action for the field, and through it 
enables one to describe the actual state of it. 

In this paper we continue the investigation, initiated 
in [3] of the reduced density matrix of a boson field in 
contact with a dispersive medium. Such a field may de- 
scribe many other situations, such as phonons in a solid. 
The resulting photonic density matrix, of a field inter- 
acting with realistic materials, is in general not in a pure 
state and should be described in terms of a mixed state 
density matrix. This situation is in stark contrast to ide- 
alized Dirichlet or Neumann boundaries, which are con- 
sistent with a field Hamiltonian. The main observation is 
that the resulting mixed state of the field is not thermal, 
which may be a consequence of the nature of the model 
we consider, which is, in essence, integrable, and thus not 
typical. 

As a measure of the field-matter entanglement, we use 



the entanglement entropy of the field. The entropy of 
radiation coupled to matter has been considered in nu- 
merous works. However, usually, the focus is on a single 
degree of freedom coupled to a bath of oscillators: For 
example, the entropy of a spin in the spin-boson model 
within the frame work of the Caldeira-Legget model [5 
was considered in [B] and the entanglement of a single 
radiation mode with an array of spins was studied in [7] 
for the Dicke model . The entanglement between spa- 
tially separated intervals of vacuum (or ground state of a 
spin chains) has been considered in [HIE]. Here we have a 
system distinct from these works, and more akin to situa- 
tions studied in macroscopic electrodynamic effects such 
as the Casimir and Van der Waals interactions between 
dielectric bodies. 

Upon considering the entanglement entropy of the ra- 
diation field in contact with a material, it becomes clear 
that the entropy suffers from a UV divergence, and is thus 
explicitly cutoff dependent. Instead, in order to look for 
universal long-range features, we consider the distance 
dependent part of the second Renyi entropy of a field 
interacting with two distinct bodies. The second Renyi 
entanglement entropy is substantially more tractable an- 
alytically the full Von-Neumann entropy, but often con- 
tains the correct scaling behavior. Recently, Renyi en- 
tanglement entropies have been computed for numerous 
systems. For example, the ability to numerically access 
1S2 has been used to probe entropy and topological en- 
tropy in recent works [TUHT2] . 

In particular, we estimate the second Renyi entropy 
and find that: 

S 2 (AUB) -S 2 (A) -S 2 {B) = 

-"U"pbVaV b ^+0(±), (1) 

where 82(A) = — logTrp^. Here p^.A is the density 
matrix of a field in contact with a body A, after the body 
degrees of freedom have been integrated out. This result 
is of interest, as it shows that the decay of entanglement 
in this model is very slow compared to typical power laws 
in quantum electrodynamics such as interaction energies 
between dielectric bodies in the Casimir-Polder regime 
and compared to the toy- model studied in 

The paper is organized as follows. After a brief in- 
troduction of the problem, we proceed to review the de- 
scription of density matrices of Gaussian states and their 
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entropy. We then use this formalism to describe the state 
of a field in material medium using the equal time corre- 
lations of the field and field momenta operators. In the 
following section, we turn to the description of the dis- 
tance dependent part of the entropy of a field interacting 
with two objects. We derive a formal expression for this 
entropy. Finally, we study the distance dependence of 
the second Renyi entropy obtaining the result ([I]). 

GAUSSIAN EFFECTIVE ACTION 

The basic principles for describing field fluctuations 
in the vicinity of heated bodies have been comprehen- 
sively explored since the early days of electrodynamics 
(see e.g. [4j [13]). Within this approach, the macroscopic 
field interaction with the material is described through 
the dielectric response function e(w,x). Following this 
logic, here we study a simplified scalar field version of 
the electromagnetic field action 

S = hj d3 -™WC( x )[ w2£ (^x) - V 2 ]^(x). (2) 

Throughout the paper we will also use the susceptibil- 
ity x, which is related to the dielectric function e by 
e(u,x) = 1 + x(w,x). 

When the permittivity is independent of uj, the action 
is local in time, and one can easily quantize the associated 
scalar action assuming the conjugate momentum 7r^ can 
be expressed in terms of (f> and doesn't depend on external 
fields. Such an action follows from the Hamiltonian H = 

J d 3 x[j^y + (V0) 2 ]. It describes, at zero temperature, 
a pure state, and as such will have no entropy. Note that 
x dependence of e does not interfere with this property: 
it just means that the field has a spatially non-uniform 
mass, but can still be describe in terms of a Hamiltonian. 

The situation is fundamentally different if e is u de- 
pendent. The non- locality of the action ^ in time sig- 
nals that our system is coupled to external degrees of 
freedom which have been integrated out, yielding a non 
trivial temporal response kernel. In such a case, the sys- 
tem cannot be in a pure state even at zero temperature, 
implying that our radiation field is entangled with the 
matter fields. 

The effective action ^ can be obtained from integrat- 
ing out other fields. To keep in mind a simple model 
for such a procedure we have in mind a bosonic field 
coupled with a matter field ip. We consider the following 
typical action: 

Spuic = 1 dt{ J d*x[cp (-0? + V 2 ) cp] + 

f B d 3 x[^ (-df - Wo 2 ) ^ + ^{d t 4>)\} (3) 

The action ^ corresponds to the form x( w ) — M '2-^2 
of the response of to a transparent, but dispersive 
medium. 



Dissipation may be introduced in a similar way, but 
requires coupling to an infinite bath of oscillators for each 
field degree of freedom, as done, e.g. in the Caldeira- 
Legget model [5]. 

When quantizing a system starting from an effective 
action such as ([2| it is important to keep in mind the fol- 
lowing subtle point: the quantization procedure cannot 
be complete without additional information on the sys- 
tem. In our case, we will need the conjugate momentum 
operators to the field and their correlations. It is impos- 
sible to obtain those from the effective action alone. The 
reason is that while the action gives us full information 
about the field correlators, it does not define uniquely the 
momentum correlations: The momentum correlations are 
extracted from the time dependence of the field correla- 
tions through equations of motion, but these depend on 
the particular way the field is coupled to the matter. 

To illustrate this point, we show that there may be 
some ambiguity on how to correctly choose those. For ex- 
ample, the terms (pdtip and —ipdt4> m a field Lagrangian, 
while classically the same, as they are related by a full 
time derivatives, yield the same effective action for <p 
upon integrating ip out but, are not quantized in the 
same way. Indeed, let us look at the following simple 
example. Consider the two Lagraniangs: 

Lx = \ (j> 2 - ^o V) + \i? ~ + a^cp (4) 

and 

Li = \ (i? - "o V) + \i> 2 - B^ 2 - aH (5) 

which differ by a total derivative: L\ — L2 = a-^{(pip). 
Let us check, that when canonically quantizing them we 
get different Hamiltonians: From the Lagrangian L\, we 
get the canonical: P$ = <f> and P^ = ip + a<fi While Li 
gives us: = 4> — a and P^ — ip. 

Clearly, the time dependent correlation functions 
(4>(t)4>(t')) computable from the action ([2]), will yield dif- 
ferent results for (tt^^tt^v)) when computed from the 
equation of motion obtained from of the full Lagrangians 
L\ and L^. 

In this paper we proceed choosing 7r^ = <p motivated by 
the usual lagrangian coupling the electromagnetic field A 
with material fields through A ■ J. 

GAUSSIAN STATES AND THEIR ENTROPY 

For a general many-body density matrix, determina- 
tion of the state requires, in principle, the knowledge of 
all matrix elements of the density matrix. These scale 
exponentially with the number of degrees of freedom and 
for interacting systems become intractable very quickly. 
However, fortunately, for states described by Gaussian 
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field theories the situation is considerably simpler. In- 
deed, by virtue of Wick's theorem, all correlation func- 
tions, and thus all matrix elements can be obtained from 
the two point functions of the field. Thus, our task is 
to use the two point functions in order to represent the 
general state of the field. To proceed, in this section wc 
review the method of calculating entropies of Gaussian 
states (see, e.g. OE]) from correlation functions. The 
material is standard, but for convenience, we choose re- 
view it here in detail. 

Consider a scalar field with n degrees of freedom cj> n 
and conjugate momenta 7r„ = —id^ n . It is convenient to 
bunch these together, defining the vector: 



(Oi,..0 2 J 



■<t>ni ^l) ■■ 7r n)- 



(6) 



Using this vector, the canonical commutation relations 
[4>m 7i"ro] = i&nm may be expressed as: 



[Oj,O k ] = iu 3 ,k 
where a is the 2n x 2n block matrix 

/„ 

-In 



(7) 



(8) 



and we set H = 1 for convenience. 

The commutation relations ([7| as well as the hermitic- 
ity of the canonical field and momenta operators are 
preserved under the group of symplectic transformations 
Sp(2n), the set of real matrices W such that WaW T = a. 

As mentioned above, we would like to use the two point 
functions of the field to determine it's state. To do this, 
let us first define the the covariance matrix: 



ljk = 2Re((0 J - (0 3 )){O k - (O fc »>. 



(9) 



The operators Oj can always be redefined so that (Oj) — 
0. In the rest of the paper we will assume our Gaussian 
field and field momenta have been correctly defined to 
have this property. 

Next, we bring 7 into the canonical Williamson normal 
form: 



\v-\V< =( I ° W 



(10) 



Where p, is a diagonal matrix jl = diag(p,i, .../i n ) an d 
W G Sp(2n) . The //, are called the "symplectic eigenval- 
ues" of 7, and they are equal to the positive eigenvalues 
of ivy. 

To see this last property, let us check that if iajip — fitp 
then /i 2 is an eigenvalue of M 2 . To see this one can 
explicitly check that uMuM — —M 2 , therefore: 

M 2 = -aMaM = -aWjW T aW-fW T = 
~W- T W T aW 1 a 1 W T = ~W- T a 1 cr 1 W T . 

We conclude that 

M 2 (W- T if>) = W- T (ia~/) 2 W T W- T ip = pL 2 {W~ T ij){ll) 



An additional, technical simplification, occurs, if we 
assume that no (<jyir) correlations are present, i.e. (4>tt) = 
0. In this case the symplectic eigenvalues are equal to 
square roots of eigenvalues ofQW, where Q and T~L are field 
and field momentum two point functions, respectively. 
To see this property, consider a 7 which is block diagonal 
of the form: 



7 = 2 



Q 
H 



(12) 



Assume that ia'yip = flip. It follows that ■y 2 tp = /J, 2 ip. 
Computing, explicitly, 



Thus 



(^) 2 ( t 



^ 2 = 4 (T QH 



(13) 



H2KSJ^H£J (14) 



and we conclude that /j, 2 is an eigenvalue of QTi (and, 
equivalently, oiHG)- 

Next, we find the Gaussian state associated with the 
covariance matrix 7 by comparing to the covariance ma- 
trix of a general quadratic Hamiltonian in a thermal 
state. 

Consider a quadratic Hamiltonian for coordinate and 
momenta cj>, tt 



H = tp Mip 



fl> = 



(15) 



where M is symmetric and real. We find the normal form 

AO 



S T MS 



and let: t/> = Sip. Then 



A 



(16) 



(17) 



In terms of the components </>i,7Ti of ip, the Hamiltonian 
breaks into independent oscillators of the form: 



H — YiiHi ; 



If, : A, ; (fa 2 + = 2Ai (afai + \) 



(18) 



where = \ S>i — iffij are canonical bosonic creation 
operators. 

Let us compute the finite temperature expectation val- 
ues of 4>i,TTi. At an inverse finite temperature /3 we have: 

'$) = h{(fli + + h) 2 ) = 

i (25,+a, + 1) = fii + § = e ^ A 1 ._ 1 + 1 (19) 
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where n, = a,i + ai is the number occupation of mode i. 
Thus: 




1 + - 

' 2 











Finally, we rotate to the original basis, so that 



(TTi4>j) (TTiTTj) 



2 S 1 







1 l S 



(20) 



Comparing ( 10 ) with ( 20 ) , and setting (3 — 1 we see that: 

Mi 1,1 
2 



2 2A, _ j 



Note that the choice of f3 is arbitrary, since one can 
always rescale simultaneously /? and A while keeping 
their product constant. Inverting the relation we find 
e 2A„« _ 1 = 2 ( /i( _ and fi na ii y: 



2Ai = log 



Mi 
Mi 



(21) 



Thus, we conclude that given a covariance matrix 7, 
and assuming this is a Gaussian state, then we can de- 
scribe the state of the system as: 



= z- 1 e -£« aA *° +<<M : 2Ai = lot 



Mi + 1 



Mi 



(22) 



Next, let us compute the entropy of this state in terms 
of the symplectic eigenvalues \ii . Since the density matrix 
factors into independent modes, we may consider first a 
single i mode. For a non interacting bosonic state defined 
as pi = -g-e -2Ai ° iCLi , the eigenvalues are labeled by the 
number of bosons occupying the mode. The occupation 
probabilities are: 



1 

1 



Si 1 



(23) 



where for n = 0, ..00 and we defined fj = e 2Ai . Let us 
compute the normalization: 



Z; = Tre 



The entropy is: 



-2A 



00 1 

\ pn 



n=0 



"6 



(24) 



= - E^O Pi,n log ft, n = - E^=0 % lQ g (£) 



logZ, - ^ E~ o<? = logZi - S^%^fe) 



-iog(i-eo-^gy 



(25) 



Substituting & = e 
of the mode is given by: 

Mi + 1 - Mi + 1 



we find that the entropy 



S, 



1 , Mi 

-log — 



1 , M* 
— log — 



(26) 



Since the modes are independent, the total entropy is 
simply: 

s = Ei K^i) ; 
/ l ( M ) = i^log^-^l gM f l (27) 

Let us check how this works for the simple harmonic 
oscillator in the ground state of a Hamiltonian H = 
J—p 2 _|_ L ^- 7 s±x 2 . We know that the state is pure and so 
should have zero entropy. Here we have a single degree of 
freedom so 7 is a 2 x 2 matrix, which is easily computed: 



7 = 2 
now we have: 



[xx) 
<pp) 



)- 


rnuj 


°) 




\ 


mu J 



muj 

■^r 



(28) 



(29) 



with eigenvalues ±1. As remarked before, the symplec- 
tic eigenvalues of 7 are the positive eigenvalues of ia"/. 
Therefore the entropy is given by: 



S = h{l) 







(30) 



More generally, we can compute the Renyi entropies, 
defined as: 



1 



1 — a 



logTrp". 



(31) 



where, in particular, Si — — Trplogp is the Von Neu- 
mann entropy. 

First computing, as before, the entropy per mode, we 
have that: 



Trpf = -—Tre 



And thus: 



-2aA i a + i a i 



^a,% - !_ a iOg 



(32) 



(33) 



Expressing S a ,i explicitly in terms of the symplectic 
eigenvalues we find: 



5<> 



io g (m a -(^) Q ) 



-l + a 



(34) 



and the final expression for the Renyi entropy S a , is ob- 
tained by summing over all the modes: 



S a — Ei ha(jM) i 

Mm) - ^gMD . 



(35) 



In particular, we will be interested in the second Renyi 
entropy, given explicitly as: 



S a = E l 1 °g(Mi)- 



(36) 
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DENSITY MATRIX OF THE FIELD IN A 
HOMOGENOUS MEDIUM 

In the next sections, we explore the state of a scalar 
field using the formalism elucidated above for Gaussian 
states. To do so we compute the elements of the covari- 
ance matrix and it's symplectic eigenvalues. 

We note that carrying the symplectic transformations 
needed to diagonalize Q implied by our action |2]), in 
general may be a hard task. However, the important 
case of a homogenous system allows us to proceed ana- 
lytically since for a translationally invariant medium the 
symplectic eigenvalues can be labeled by momentum. 

In the case of a homogenous medium, the action ([2| 
becomes: 

S= lrJ (^3 d ^(k)[^,k)+k 2 ]^(k). (37) 

As explained above, the density matrix of the field can 
thus be written as: 



Pfield 



(38) 



where is a bosonic creation operator labeled by mo- 
mentum k. 

For the models we consider here, the system doesn't 
have (</>7r) correlations, thus we proceed by computing 
the field two point functions Q and field momentum cor- 
relations %. As explained above, the desired symplectic 
eigenvalues /ifc are square roots of the eigenvalues of the 
matrix QT~i. In a homogenous space, both Q and T-L are 
diagonal in momentum space, and we are left with the 
task of computing the momentum space correlators. 

We thus have: 

Explicitly, to compute the correlation functions, we 
will use that for a single harmonic oscillator, governed 
by an action 



do; 
2^ 



AH0( w )0», 



we have: 



(0(4)0(0)) 



50 he~ iujt dw 
o A{iui) 2ir 



(39) 



(40) 



Therefore, taking account that the momentum is a good 
quantum number we can write: /ifc = 27r~ 1 (gk^k) 1 ^ 2 , 
where the action ^ allows us to compute, per k mode: 



f°° 1 
9k = (M-u)=J d--^— 



(k, iw) + k 2 ' 



(41) 



Similarly, using time point splitting and the equation of 
motion ir — 

h k = (7r$) k = limt/_rt d t d v f °° dej ^ffj^ = 

r ^ Mk i )+fc . (fc 2 +x(*H)). (42) 



As usual, quantizing the fcs according to k — 
-^(ni, ..rid) where n, are integers, and L — >• oo is the 
linear size of space, we have that 

E ni ,. nj h (P^(n u ..n d )) = LdL ~ d I2n u ..n d h iP^( ni ,..n d )) 



(43) 



In other words, the field entropy per unit volume can be 
written as: 



S field oc J d d kh(fik) 



(44) 



Let us check the resulting symplectic eigenvalues \ik 
for a simple case: that of a free Gaussian field. In this 
case we have: 



'a 1 _ w 
5fc < /ree ~/ dUJ lo^¥~2{k\ 



and, 



r°° [ 

hkjree = I duj 



I 

JO 



k 2 = ^' fc| 



2 + fc 2 2 



(45) 



(46) 



we immediately get that the symplectic eigenvalues are 
all 1: 



— 2?T 1 (^k,free^k,free) 1 ^ 2 — 1- 



(47) 



Thus, we immediately get Sfi e id = since in the entropy 
formula (27), h((j,)\n=i = 0: the action describing a free 



field does not carry any entropy. 

The last calculation may seem as a rather convoluted 
way of reaching the simple conclusion that the free field 
is in a pure state. However, the calculation serves as 
an important check for us that the present treatment is 
consistent. It is also clear how it can be adapted to more 
complicated situations, as studied in the next sections. 

Next, we consider the state of the field choosing a 
typical dielectric function e to use in ( 41 1 and (42 1. 



Concretely, we use a typical susceptibiliy of the form 

, ,2 



(a 



-vy p u>) 



Using the residue theorem, the integrals (41), (42 1 can 



be expressed in terms of the roots of the fourth order 
polynomial appearing in the denominator of the inte- 
grands. However the expression is rather cumbersome, 
but can be easily used to numerically study the symplec- 
tic eigenvalues of 7 as function k. In Fig [T] we exhibit the 
k dependence of the effective energies Ek = 2A& appear- 
ing in ( [38] ) for various values of the coupling strength of 
the field to the medium as represented by w p . 

FigJT] clearly shows, that k dependence is very differ- 
ent from the typical linear dispersion of photon energies, 
showing that the state is not a thermal state of photons. 
It is also clear that the low momenta and high momenta 
asymptotics of the Ek are quite different. In the following 
sections we proceed to compute the symplectic eigenval- 
ues of the field covariance by asymptotic analysis of the 



integrals (41), (42) in the low and high k limits. 



G 



Effective energies faff (id = 20 t y = 1 




— Mp = a 



Carrying through the integrals we find again two cases: 
1) 7 P 2 > 4ll>o 2 - Then a± are real and negative and 



III 



log[fc]7 p log [-a_] a 2 



log [-a + ] a 2 . 



=(53) 



2) 7 P 2 < 4a>o 2 - a± are complex conjugates, we find: 



200 4CO 600 BOO 1000 



TTT — 7p log[fc] 2f p log[u ] 
111 _ fe 2 2fe 2 

(7p-2w 2 )arg ^ p -iy^wo 2 -~t p ' 



/c 2 v /4cd 2 -7p 2 



FIG. 1: Effective energies for various values of the coupling 3) g — In this case" 



High k asymptotic behavior of symplectic 
eigenvalues 

To compute the cutoff dependence of the entropy for 
X&i we consider the large k asymptotics of gu,hk- Using 



(41| at large k we rewrite g k as: 

l 



PCX 

9k = Jo 



w 2 + k 2 +UJ 2 {1 _ B) 



dw . 



(48) 



where B = ?* w+ "° Noting that ulB < ujI, we 
can expand in series in £>, obtaining, to lowest order in 
B/(u 2 + k 2 + u 2 p ): 

o, - f°° dw( 1 Vu? x - ls^±4 ) 

higher orders 



2 7^+fc 2 " 1 ~P 
1 IpW+ujI 



will (49) 
let us com- 



With II = f_ dwUjl 7-1 T" r-f : — r 

J0 P (m 2 +fc 2 +t<; 2 ) 2 u) 2 +7 p i«+tiJ5 

pute II. 

Writing: w 2 + 7 p cj + = (cu — a_) (w — a+) , where 

a± = 7p=l=v ^ 7p ~ 4a ' , there are two cases: 

1) 7 2 > 4w 2 . In this case a± 
are real and negative and carrying the integral we get 



tt _ 7p log k _ 7j_ . 
11 _ fe 4 2fc 4 + 

- logf-a^la^+logj-a+jq^ 
fc 4 (-a_+a + ) 



(50) 



2) 7 2 < 4wg. In this case a± are complex conjugates 
and we find 

tt _ 7p log k _ 7p+7p log^o , 
11 ~~ fc 4 2fc 4 + 

fcvfef ~ 2 ^ ar ^ - i v 4 "° 1 (51) 



In particular, if 7 P = we have: iT 
we write: 



Similarly, 



/ife = (fc 2 
where to lowest order: 



w p).9fe -w 2 III 



(52) 



III 



+ IpW + uJq w 2 + k 2 + i 



-dw 



III = 



2fc 2 



(54) 



(55) 



Computing to lowest order in ui p the corrections II, III, 
we find for large k, small lo p and 7 P 7^ 



tt 2 w 2 7 p 7T log k 
9khk " T + 2P ■ 



(56) 



We remark that if 7 P = we get a slightly modified result: 

.9fc/i/c - J^ — jp— • 

We conclude that the symplectic eigenvalues /ik, to 
lowest order in oj p , behave as: 



w 2 7 p log k 



+ .. 



(57) 



for k » 1. 



Low fc behavior of symplectic eigenvalues 

Studying the low k behavior requires slightly more care 
in the expansion. Here we find that the leading behavior 
for the field correlators at low k is given by 



9k 



2k J I 



(58) 



To do so, we start by rescaling w in the integral (41 1: 

1 



roc 

9k = Jo 



-dw — 



w 2 -\-k 2 -\-- 



1 rOO 
k Jo 



l+a; 2 +:E 2 



dx. 



(59) 



We continue to split the integral into two intervals 
(0, fc -1 ) and (fc -1 ,oo), and separately approximating 
each term: 



k JO 



l/k 



-dx 



l+x 2 +x 2 



Hi- 



kx\kx-\- g ) \ 
fc 2 ai 2 +7pfcrr+^ 2 > 



1 r 
fc Ji 



1 f l/fc 
fe Jo 



fcx(fcx + 7p ) 
t 2 - 2 +Tpfe=: + ^' 



l+x 2 +x 



2 ^( 1 --, 

^ + \ Ii/k TTx 2 ^ 



(60) 
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for small k we get: 

1 

k 



-dx -» 1 



l/k 



(61) 



and: 



k Jo 



l/fc 



z?-d x 

l +x 2 +x 2P 



2 k 

a similar treatment shows that: 



1+- 



1+- 



(62) 



o(fc). 



(63) 



Finally, we find that the symplectic eigenvalues diverge 
at k — > since 



2\2 



1/4 



2fc I K 2 +^)( 4w 2_ 7 2 +4w ^ 



(64) 



diverges as -^== for small fc. The entropy function 



for small fcs diverges as /i(fc) ~ ~log{k). However, in- 
tegration d D k is finite for any dimension and so there 
are no infra-red divergences in the total entropy per unit 
volume. 



Interpretation of the field reduced density matrix 

Having evaluated the symplectic eigenvalues of the 
field covariance matrix we can write the density matrix 



of the field as (22| 



= Z- 1 e-^ 2Aka+kak : 2A fc =loe 



Xg k h k + 1 



i9kh k - 1 



(65) 



As mentioned in the introduction, the first natural ques- 
tion to ask is: is the state of the field p thermal? This is 
certainly a natural possibility suggested by considering 
the bulk material as a thermal bath for the field. In fact, 
one may always write the state p, formally, as thermal, 



p = Z- l exp(~pH eff ), 



(66) 



for some suitable operator H e ff. However, to actually 
interpret the state as thermal, we would like H e ff to 
represent a reasonable, local, physical Hamiltonian. For 
the Gaussian states we can write H e ff explicitly as: 



H e ff = J dxdx'u(x — x')4flj 



(67) 



where 



2/T 1 / dkA k e 



ik-(x— x' 



(68) 



The locality properties of this Hamiltonian thus depend 
on the locality of the Fourier transform. We find that 
generically the Fourier transformed u gives us a non-local 

H eff- 

Alternatively, we may interpret 2Ak as 2Ak = f3 k uj k 
where w k are the photon energies in a free homogenous 
space, without the interaction with the material, i.e. 



(3 k uj k = 2A fc ~ log 



^29kh k + 1 



gkh k - 1 



(69) 



As is clearly shown in fig[T] if we take w k cx |fc| we find 
that different momenta k feel different effective temper- 
atures. 

In the infra-red limit expanding this expression for 
small k using we find the leading, small k behavior: 



f3 k u> k — v\l k 



2, 2\2 



K+^)(4^_ 7 2 +4 ^) 



-1/4 



(70) 



If we assume uj k cx |fc| we conclude that the effective 
temperature of this k mode is: 



T k — — oc \fk as 
Pk 



0. 



(71) 



What is the energy and number of occupied soft modes 
per unit volume up to a given k m ? We show that these 
are proportional to k^t 1 ^ 2 and k? n 1 , respectively, are 
finite and small. Indeed, since the occupation of a mode 
k is given by: 



1 



n k 



(72) 



The expected occupation number for modes with k < 
fcmin is thus given by: 



A^kmin = / d kn k oc / dk 

'fe<fe min Jo 



d— 1 



k d-l/2 



lf3) 



v\fk v [d — _ 
The energy, assuming u k oc k, is: 

E (k < k miD ) ~ J k<kmin d d kn k k cx 

JO aK uVk u{d+±) [ ' ' 

Finally, the number variance of modes up to k m is com- 
puted as: 

SN 2 (k < fc min ) = J k<kmin d d kn k (1 + n k ) oc 



f k -dkk 



d-1 J_ 
<j, 2 k 




Note that for d = 1 we find an infra-red divergence: 
(SN 2 ) = —\ogem, where em is an infra-red cutoff, in- 
versely proportional to the system size. 



Cutoff dependence of the entropy 

fn this section, we compute the total field entropy, and 
show that it suffers from a UV divergence. This result is 
not surprising, since, in principle, entanglement can get 
contributions from all momentum scales. 



Using (44), the full quantum von-Neumann entropy 



of the field per unit volume is given by the momentum 



integrals of h{\ ^gkhk)- The large momentum behavior 



of the integrand is obtained from (561 to be 



^ p "fp l°g fe(l+3 log fc— log 



"ply '°sfc 



(76) 



for a transparent medium (j p << f) we consider addi- 
tional terms and find: 

h( r^T^ ^o(l + 31ogfc-log^) 
h{\ -,g k h k ) — ; (77) 



using the integrals J lo ^ 



4fc 3 

log[fc] 2 
2 



and /33jffiL = 



log[fc] J 



we can estimate the integration of the expressions 



( 76 1 , ( 77 1 over k in 3d, to find the entropy per unit vol- 



ume 



w| 7p log(A) 3 
S field oc I w^olog(A) 2 7p = , (78) 
plasma model 



where A is a high momentum (UV) cutoff. 



Note that in obtaining (78) is obtained to lowest or- 
der in ujp. However, this approximation is justified for 
our purposes since it becomes exact at the large k limit 
which we are studying, fndeed, numerically, the approx- 



imations used in (78) actually recover the correct cutoff 
dependence even for large values of uj p , since the dielec- 
tric response decays at large k values. 

It is interesting to observe the special place of the "pure 
plasma" limit response function. We can easily under- 



stand the result (78) as follows: Substituting the plasma 

permittivity limit form \ = — in the action ([2]), in 
a homogenous space, we see that the role of x is simi- 
lar to producing a mass term for <j). Thus, the resulting 
action is consistent with a Hermitian field Hamiltonian, 
and as such, at zero temperature, to a pure state. One 
can also understand the vanishing of entropy for pure 
plasma as follows: Consider a slab of material, with a 
pure plasma form for the dielectric response. Such a ma- 
terial will have no losses: A pure plasma system will 
be completely transparent to radiation above the plasma 
frequency, and completely reflecting at lower frequencies. 



S R (A,B) = S 




FIG. 2: A Casimir entanglement entropy. 



The field will have a finite mass inside the region occupied 
by the plasma and will not generate entropy. It is inter- 
esting that the use of a pure plasma in the computations 
of Casimir energy and Entropy has been at the heart of 
a recent debate [T6lU9j . We do not make any claims re- 
garding the debate, but note that the distinction between 
the Casimir entropy in the two models is manifested in 
the full quantum entropy computed herein. 

We remark that it is possible to incorporate the high 
momentum cutoff more naturally by taking into account 
the spatial dispersion in e. While different functions may 
have different asymptotic properties, we expect to get 
similar behavior. Indeed, if we take, for concreteness 
sake, the simplest extension of the previous treatment, 
we can use (see e.g. [2D]). 



e(iuj, k) = £ ( 1 + 



A \ l 



f 

"fuj 4 



(79) 



where the expression is valid for small k <C n/a where 
a is the interatomic distance. One can easily check that 
this form doesn't change much the cutoff dependence. 



DISTANCE DEPENDENT ENTROPIES 

Since the entropy Sfi e id is UV divergent, it is natural 
to ask, in analogy with the Casimir effect, what is the 
distance dependence of the entropy of interaction with 
two distinct bodies A and B, and is the distance depen- 
dent part of it is UV finite. To answer such questions, 
the "Casimir EE" was defined in [5] as 

S R . a ( A, B) = S a ( A U B) - S a ( A) - S a (B) , (80) 

where S a (M) is the a— Renyi entropy of a field interact- 
ing with a body M described by e = Sq + \m, where 
the susceptibility xm = everywhere outside M. The 
situation is illustrated in Fig(2] 

The following few remarks regarding Sr are important 
to note: 
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1. Sr is distinct from the more familiar Casimir 
entropy. Casimir entropy is denned as Sc(A, B) = 
— limr->o 9tFc, where Fc is the Casimir free energy, 
obtained by subtracting all distance independent terms 
from the free energy of the EM field in the presence of 
bodies A, B or boundary conditions. 

2. Sr is different from the "relative entropy" of prob- 
ability theory, as we are comparing different systems, 
and not merely different statistical information about the 
same system. 

3. It is important to note that, while the sub-additivity 
of von-Neumann entropy [21] shows that the entropy S 
of a field is always larger than the total entropy of the 
combined system of field matter, this is no longer evident 
once subtractions are taking place in order to define Sr. 

The relevance of Casimir entropy Sc to understand- 
ing thermal corrections of the Lifshitz formula has been 
pointed out in many papers (see e.g. [E3 [52 [23] ) , where 
it was noticed that as T — > 0, Sc may not go to zero 
when using the Drude model, as would be expected by 
the Nernts theorem. It is quite interesting to note that 
while the Casimir entanglement entropy Sr is distinct 
from Sc , a similar behavior is observed in ( 78 ) . Sc has 



a clear thermodynamic meaning, especially at high tem- 
peratures, where the Casimir force is entirely entropic 
[2"1| . Indeed, at high temperatures we expect Sr = Sc, 
as most of the field entropy will be thermal (Technically, 
the relevant Green's function gets its major contribution 
from the oj = Matsubara pole). 



FORMULAS FOR THE DISTANCE DEPENDENT 
VON NEUMANN AND RENYI ENTROPIES 

In this section we describe the derivation of an abstract 
formula presented in [3] for the Casimir entanglement 
entropy. In addition, we generalize the formula to also 
describe arbitrary Renyi entropies. It was shown in that 



3 



a/2+ioo , '° s [^tJ 

1 /9_!~ UJ - 



Trlog 



(K A K B + K Alm -Ka-K b ) j^- )(81 



l-K, 



where Ka — a ,_ 1 j" 4+ ,^ (Ta — 1/4), and similar expres- 
sions hold for K B ,K AUB . Here T = QaH a , where Ga 
and Ha are field and field momentum two point func- 
tions, respectively, in the presence of body A. 



In this section we repeat the derivation of (81) and 



generalize the expression to general Renyi entropies: 



S 



R,a 



1 ,1/2+ioo 
it Jl/2-ioo ^ 



((2^+1)° 1 -(2^-l)° l )c 
/2-ioo Ul * ((2 v ^+l)°-(2 v / i-l) Q )v / S(«-l) 

(K A K B + K AlJB -Ka-K b ) T ± r 



The expression (81 1, is similar to the TGTG formulas in 
it's form: an integral over the TrLog of a combination of 
Green's functions. However, it differs from such formulas 
in three major aspects: 

1) The integration variable x is not a frequency vari- 
able, but rather an auxiliary spectral variable, 

2) The presence of the term Kaub does not allow for 
full separation into local object properties and free prop- 
agators. 

3) The non-analyticity of the integrand at x = 1/4. All 
of these make the formula harder to use than the TGTG 
formulas. Nevertheless, it can be used as a starting point 
for various expansions when the bodies are weakly en- 
tangled with the field, so that \\T A - 1/4|| < 1. 

The derivation bellow of eq. (811 follows and adapts 
the approach of [55] to the scattering formalism for 
Casimir energies. 

To compute the entropy, we first need to find an ex- 
pression for the density of symplectic eigenvalues of the 
covariance matrix 7. In the absence of (</>7r) correlations, 
these symplectic eigenvalues (up to a factor 2 in the def- 
inition of the covariance matrix eq. ([9])) are related to 
the square roots of eigenvalues of T. 

Thus, we first find a convenient representation to the 
density of density of states of T = QH. Note that T is 
not Hermitian, however, since Q and % are positive Her- 
mitian matrices, it has the same spectrum as Q^-^HQ 1 / 2 , 
and one may safely use the formulas below. Consider the 
representation: 



5(E - x) = ilm F * = 
±Imd E Log(E -x + iO) 



(83) 



Then we have: relative density of states of a Hermitian 
operator X as: 

5 PT (E) = pr{E)-p Ta {E) = 
ilm<9 s [TrLog(£ — T + iO) — TrLog (E - T + i0)\ = 
iTmdB TrLog ((E-T + iO)(E -T + iO)- 1 ) 



= ilm9 B TrLog ( 1 + (T - T) 



s-r +io 



(84) 



Denoting: Go = B -r a +is anc ^ ^ = ~ To), we can 
write (84) as 



Spr(E) = -Im<9 B TrLog (1 - G D) 
To compute 

Sr(A, B) = S(A li B) — S(A) — S(B), (85) 
we need the relative densities: 



i-K A 



Trlog 

(82) 



Sp R = Sp rAUB - S PTa - 6p rB = 
ilm5 B [Tr log (1 - G D A ub) ~ Tr log (1 - G D A ) 

-Trlog(l-G D B )]. (86) 
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Here Da is computed using Fa = GaW-a, where the cor- 
relations Ga,Ha are computed for the field in the pres- 
ence of body A, and similarly for Daub, Db- We com- 
bine the terms using the following identity: 

1-G D A i 1 ~ G o£>AUb) ±_g qDb = 

l= ^((l-G D A )(l-G D B )- 
GoDaGqDb - Go (-Daub — D A - D B )) 1 _ G 1 ||Jij = 
1 - i-g d a (G DaG D b + D AVB - D A - D B ) 1 _ g 1 oDe 



Defining: 

Ka 
we obtain 



GnD 



1 



0-^A 



E-To + is 



77 (T, 



Sp(E) R = ±Imd £ TrLog(l - 
^ (K A K B + K AVB -K A - K b ) ^ 



We now use this spectral density to compute the en- 
tropy using the eigenvalues of the covariance matrix, eq. 



(271, 



^logV 



The lowest eigenvalue QT-L is always larger or equal to 
1/2, by the uncertainty relations. Thus we may write 



S 



1/2 



h (2y/x) 5p(x)ndx 



(87) 



Integrating by parts and moving the contour integration 
to the imaginary axis, we write this expression as: 



S 



IT Jl 



1/2+ioo 



da;- 



TrLog 



1/2-ioo 

(* - T=fe ( X ^^b + ^aub -K A - K b ) j^j 



which is eq. (11). 

We can also extend this formula to cover Renyi en- 



tropies (31). Using the formula for the density of states 



(861 together with the expression for the Renyi entropy 



in terms of the symplectic eigenvalues (34), we find 



1/2+ioo , 
L/2-ioo ' 

1 {K A K B + K hVB - K A - K B )j^K-)dE (88) 



■ltij£r{d E h a WE))Tn*g{i 



1-K A 



yielding ( |82[ ). In particular, the second Renyi entropy is 
given by: 



S 



R,2 



IT Jl/2-ioo UX 5-r- 11 1U &^ 



l-K, 



(KaKb + Kaub-Ka-Kb)^-). (89) 



To relate this form to TGTG formulas, we recall that in 
such formulas, the relative density of states of the elec- 
tromagnetic field interacting with for two bodies A, B, 
through dielectric susceptibilities Xi(w),xs(w) at fre- 
quency u} is expressed (before Wick rotation) as: 



1 



Imd u Tr log(l - TaG q TbG ) 



(90) 



where Ta — ^> 2 Xa i-Gquj 2 xa are ^ e Lippman-Schwinger 
operators of the problem, and Go are free propagators. 
Choosing Ka 1 J L Ka to play the role of TaGq and 



Kg 1 _ Kb the role of TbGq, The density of states ([87]) 
may be written as: 

Sp(E) R = 
ilma £ TrLog(l - T A G T B G Q 
1 (K AU B-K A -K B ) 1 



X-K A 



l- Kb 



We observe, however, the appearance of an additional 
(Kaub — Ka — Kb) term, which is not-separable into a 
product of correlators of the separate bodies. 



REPRESENTATION OF THE CORRELATION 
FUNCTIONS IN TERMS OF 
LIPPMANN-SCHWINGER OPERATORS 

In the computation of the correlation functions bellow, 
we will use extensively the representation of correlation 
functions in terms of Lippmann-Schwinger operators. For 
a body A, we define the Lippmann-Schwinger operator 
Ta at imaginary frequency iuj by 



T A = uj xa 



1 



1 +ff0W 2 XA(«w) 



(91) 



where go are free propagators. In particular, for a scalar 
field we take go = _ A 1 f ^ 2 ■ 

The Lippmann-Schwinger operator Ta is related to 
the green's function gA , 



(x\g A (w)\x') = (x\ 



1 



-A + w 2 £a(«w) 
by the operator equation: 

g A (uj) = g (oj) + go{uj)T A go(u) 



\x') 



(92) 



(93) 



The eq. (41) for the (<jxj)) correlation functions written 
in a general basis (i.e. without assuming translational 
invariance) is then: 



Q a = J °° dug a = J °° dw (g + g T A g ) 
So + Jq 00 dujg T A go 



(94) 



Similarly, as in ( 42 ) , the field momenta correlation func- 



tions are encoded by 
Ha = Wo 



&uju 2 goT A go- 



(95) 
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(in the expressions above, and what follows we omit the 
lo dependence in go). We note that in the continuum, 
Go, Ho are diagonal in momentum, with matrix elements 
given by ^pry S(k — k') and ^^-S(k — k') respectively (as 



obtained in eq. (|45|) , (|46| ) . 

Let us recall some general properties of T A on the imag- 
inary frequency axis |25) . Below we will repeatedly use 
that, as a consequence of the Kramers-Kronig relations, 
combined with the assumption of equilibrium we have 
XA{iu) > 0. In it is known that xa{iw) is real and de- 
caying as uj~ 2 as w — > 00. 

We have the following properties: 

1. T a i/j = for any tp for which vanishes on A. This 
is established by rewriting Ta as 

1 



T A = 



Xa 



1 + ^/uj^xa9oV^Xa 



\[^XA (96) 



2. Ta is a positive operator, i.e. (%P\Ta\iP) > for 
any tp in the Hilbert space Ta acts on (square integrable 
function supported on the region A). 

And, 

3. 



T a <oj 2 (£aM-1) 



(97) 



as operators, i.e. for any vector ip, (ip\T A \ip) < 
(^ 2 (e A M-l)|V>. 

In the next section we use these properties in our anal- 
ysis of the distance dependence of the entropy Sr. 



LARGE DISTANCE EXPANSION OF THE 
RENYI ENTROPY S 2 ,r 

In this section, we study the behavior of the relative 
entanglement at large distances. Consider two bodies 

with a dielectric function ea{ iu) — l + ^ 2 ^ 2 (and similar 

expression for body B), and volumes V A ,V B . 

Our main result is that at large separation the Renyi 
entropy of the field is the sum of the separate body en- 
tropies, with the correction decaying as: 



2 2 27r 
S 2 ,r = ~u pA u pB V A V B 2 4 



(98) 



Note that the power law differs from the typical power 
law of R~ 7 appearing in the Casimir-Polder interaction. 
To find this result we first write the 2-Renyi entropy 



(36) as: 



S 2 = Slog M = -Trlog(l + 5r) 



(99) 



Thus the relative S 2j r is given by: 
S 2 ,r = |Trlog (1 + ST B ) - 1 (1 + 6T AVB ) (1 + 5T A ) - 1 



i+sr. 



iTrlog 



(<5r A uB - ST a - ST B - 5T A ST B ) 



1 



(100) 



and in particular, in the long distance expansion, we ex- 
pect: 



\\sr 



AUB 



-ST A -ST B -ST A ST T 







(101) 



and we can approximate: 

S 2 , R ~ \tt (sr Al]B ~ 8Y A - sr B - 5r A dr B ) (102) 



In calculations of (102), we have several different kinds 



of terms. We concentrate on the so called "dilute limit 
" where it is assumed w 2 x(jcj) <C 1. In this case we can 
use the approximation Ta ~ u 2 XA(iw). 

To lowest order in \A, Xb, we have the following terms: 

S 2 m ~ 5Tr("H (^AUB - SQa - 8Gb) + 
(SHavb - SH A - 5U B )Go - 
HoSGaHoSGb - SHaGoSHbGo)) (103) 

It turns out that the leading contribution is obtained 
from: 

Tr% SgAH Sg B + 5H A GoSn B Go = O (iT 4 ) (104) 
The calculation goes as follows. Explicitly, using 



(451,(46) ,(94) and (95 1 



TrHo6G A no5GB + $U a Qo5UbGq = O (R- 4 ) = 
J dujduj 1 ' d?k(k\Ho9oT A goHo9oTBgo + 
u 2 ui' 2 Go9oTAgoGogoTBgo\k) = 
J dujdu' d 3 kd 3 q k2 ^ 2 T A (k,q,u) g2 ^ 2 x 

i^lMg.fcV)^ + |fc||«l) (105) 

Now consider the effect of shifting the object B by a 
vector Rh to Br — {x : x — Rh € B} . The Lippmann- 
Schwinger operator Tb r associated with the shifted body, 
written in momentum representation, is 

T BR {q,k,uj') = J dxdye i ( kx - q rt T B (x + Rn,y + Rh,u') 
J dxdye 4 ( fe - x -9-^- i ( fc -«)- J?ii T B (a;,y,a;') = 
e-^ k - q > m T B (q,k,ij') 

We therefore write for the shifted position: 

K(R) = 

JdujdJ d 3 kd 3 q w ^T A (k,q,uj)-^ j2 x 
e^^T B (q,k,u>')^(^ + \k\\ q \) (106) 

We now rescale all momenta and frequencies appearing 
in the integral by: k = Rk. 

K = d 3 kd 3 a &^T A (i f , |) x 



1 e-'(fe-g-)- 
q 2 -\-u) 2 q 2 -\-tj' 2 



q k u)' 
R> R> R J fc2 + . 



fc 151 
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Note that for R -» oo we have Tb |, 4, | 



(I- 



Note that F = at i; -> oo at g -> oo as well as 



since in this limit x ~^ an d 
T B can be approximated as in the dilute limit. For 
concreteness, let us take 



F(0,0) = 1 and £2(0,0) 



Thus we have: 



1 x e A 
otherwise 



V A V B y£ pB I r °° dkdq(d k d q F) cos(fci? Wo ) cos(qBu ) 



A (z) = 



+2 / °°dk(d fc F(/c,0))cos(/ci?u;o) 



(116) 



(107) 



We can now carry out the frequency integrals yielding: 



The remaining integrals in (116) can straightforwardly 
be shown to decay as Rr 2 (using more integrations by 
parts), showing that: 



K 



1 k\ e -i(k-q)- 



ir 2 (2k 2 q 2 +2kq{k+q)Ru) + (k+q) 2 R 2 u:o 2 ) 
Akqik+qYik+RuJo^iq+RuJa) 2 

For R —¥ oo we can use the approximation 

= / dxe^'^-?) ~ V A 

J A 



K = V A V B - 



u> 2 R* 



+ 0(4r), 



(117) 



(108) establishing the asymptotics (98 1. 

The complete analysis we check that the additional 
terms in (103) give a sub-leading correction to (98 1. In- 



fo q 
R' R 



deed, the term 



(109) 



where Va is the volume of body A. We therefore have: 
K-VaVb 7 ^^ J d Z kd 3 qx 

(2k 2 q 2 +2kq(k+q)Ruj + (k+q) 2 R 2 u> 2 ) _ l(k _ q) . fl . . 
4kq(k+q) 2 (k+RuJo) 2 (q+Ruio) 2 C ^ W> 

where we also used that \go,Ho] = 0. Using (46) for Ho, 
We can now carry out the 3d angular integrals in the we nave 



TrH (<5<?AUB ^ SQa - 5Q B ) 
= TrHo J dujgo (T AVB - T A - T B ) g (118) 

gives us a R~ 6 decay. The analysis goes as follows. Using 
cyclicity of the trace this expression is the same as: 

J^duTrpAVB-TA-TB) g H g (119) 



standard way, writing: q ■ n = cos 8\q\, we get: 

(2k 2 q 2 +2kq(k+q)Rl» + (k+q) 2 R 2 u> 2 ) ■ /, s . / x , . x 

( k+q )z( k+RujQ )2( q+Ruo) 2 sin|Kj sin v q; {ILL) 

At this point, it is convenient to rescale back the mo- 
menta, writing: 

K = u? pA u? pB V a Vb ^ / °° dkdq x 

(2k 2 q 2 +2kq(k+q) + (k+q) 



(x\g (u)n g (uj)\y) 

j d 3 k |,| 1 i(x-y)-k _ 

J (27T)3 l K l(fc2 +tJ 2 ) 2e - 

A f n °° dk f\ dX n2 f 2y2 e^~y\ Xk = 

47T" ! JO J— 1 (fc 2 +W 2 ) 

2*4-1/1 Jo" dk (k4l 2 ) 2 Sin (' X ~ ^l fc ) = 



(fe+g)2(fe+l)2( 9+ l) 



sin(fci?w ) sin(gi?cjo) (112) 



To analyze this integral we do a couple of integrations by 
parts according to: 

ff Q L dxdyF(x,y)d x d y G(x,y) = (F(x,y)G(x,y)\^ ) |£ =0 
J L dj(d y F(x,y))G(x,y)\^ - 
tfdx((d x F(x,y))G{x,y)\i =0 ) + 
J Q L dx J Q L dy (d x d y F(x, y)) G(x, y) (113) 

In our case we take: 

(2k 2 q 2 + 2kq(k +q) + {k + qf) 
~ (k + q) 2 (k + l) 2 (q+l) 2 



In the dilute approximation: 

/ °° duTr (Taub -T a -T b ) goH go ~ 
Jo° dw Iaxb dxd y w4 XA(«w).goAB(a;, y)xs{ioj) x 
(x|g (w)'Ho5o(w)|y) = 



(120) 



(« 2 +l) 



* sin 



x-y\ 2 



2tt) 3 \x-y\ 6 Jo duJ XA (i \ x - y \ J XB (i \ x 



-v\ 



Io° du (u2 U +L) 2 sin MV 



3„-ui 



' R— j-oo 



(114) 



Lx B dxdy^™^ /rdwx 

2 , ,2 



and 



5 



(27r) 3 W(! JAxB UXU -y |x-y| 6 



G 



■ cos(fc Rujq) cos(qRuJo). 



(115) 



giving us a contribution of 
sub-leading to the 0(R~ 4 ) contribution (104). 



5{2ir) 3 R? 



1240.R 6 ' 



(121) 
which is 
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DISCUSSION 

In this paper we continued the investigation initiated 
in [3] of the state of a field interacting with a dispersive 
medium within the Gaussian model. We showed that 
the state cannot be considered as thermal, but rather 
as a state where photons have a k dependent effective 
temperature. We found that the field is described by a 
density matrix whose Von-Neumann entropy diverges as 
described by eq. ( [78| . 

In addition to supplying details on some of the calcu- 
lations carried out in [5] we present several new results: 
Namely, formulae for the distance dependent Renyi en- 
tropy , as well as the distance dependence of S2,r- We 
find that the decay in 3d is proportional to S2,r °c i?~ 4 . 

This result is curious, in that it seems at odds with 
the scaling Sr ~ i?~ 6 of the entropy for parallel plates 
per unit area found in [3J. Indeed, using the asymptotic 
result (98 1 we can approximate the distance dependent 



part of the entropy of two plates per unit area by a pair- 
summation as: 



S 2 ,r oc J dxdy 
2tt /"rdr 



(R 2 +x 2 +y 2 ) 



77^ 



(R 2 +r 2 ) 2 



(122) 



giving us a much slower decay compared to that found in 
[3J . It must be noted that the exactly solvable toy model 
considered in [3] is unusual in that it has a ip<j> coupling 
rather than the <fnp coupling considered here (See the dis- 
cussion in section II), but it is not clear if this difference 
is the source of the different scaling behavior. Also, the 
calculation in [3] was done for the full Von-Neumann en- 
tropy rather than S^,r considered here. More work is 
needed to understand the difference between these re- 
sults. 

We expect that much more insight into the mixed state 
of the electromagnetic field may be gained using numeri- 
cal means to study how other factors, such as geometries 
and vector properties, affect the Casimir entanglement 
entropy and entanglement spectrum. 
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